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A self-consistent theory is proposed for the general problem of interacting undulating fluid mem¬ 
branes subject to the constraint that they do not interpenetrate. We implement the steric constraint 
via an exact functional integral representation, and through the use of a saddle-point approxima¬ 
tion transform it into a novel effective steric potential. The steric potential is found to consist of 
two contributions: one generated by zero mode fluctuations of the membranes, and the other by 
thermal bending fluctuations. For membranes of cross-sectional area S, we find that the bending 
fluctuation part scales with the inter-membrane separation d as aP 2 for d <C vS, but crosses over 
to d~ 4 scaling for d y/S, whereas the zero mode part of the steric potential always scales as d~ 2 . 
For membranes interacting exclusively via the steric potential, we obtain closed-form expressions 
for the effective interaction potential and for the rms undulation amplitude a, which becomes small 
at low temperatures T and/or large bending stiffnesses k. Moreover, a scales as d for d <C \fS, 
but saturates at ^/k^TS/ k for d VlS. In addition, using variational Gaussian theory, we apply 
our self-consistent treatment to study inter-membrane interactions subject to three different types 
of potential: (i) the Moreira-Netz potential for a pair of strongly charged membranes with an inter¬ 
vening solution of multivalent counterions, (ii) an attractive square well, (iii) the Morse potential, 
and (iv) a combination of hydration and van der Waals interactions. 

PACS numbers: 82.70.Dd, 83.80.Hj, 82.45.Gj, 52.25.Kn 


I. INTRODUCTION 

Lipid membranes are two-dimensional fluids lacking in¬ 
plane shear elasticity, and the only elastic penalty they 
experience comes from changes in curvature [T|. Thus, 
they undergo vigorous thermally activated undulations, 
and a pair of such membranes brought into proximity 
will experience an osmotic pressure generated by the re¬ 
duction of phase space available for the membranes to 
fluctuate. The problem of determining this osmotic pres¬ 
sure is a classic one, going back to Helfrich [2|, and over 
the past four decades different methods have been pro¬ 
posed to treat it, employing various degrees of heuristic 
argument, field theoretic tools, and functional renormal¬ 
ization group (FRG) techniques (see, e.g., Refs. [3lfl3]). 
A major difficulty that has dogged research in this area 
is associated with the proper treatment of the steric con¬ 
straint, i.e., the constraint that a pair of membranes may 
not penetrate each other. One response was to replace 
the problem of implementing the steric constraint at the 
level of the partition function by an Ansatz that the root 
mean square (rms) fluctuation amplitude (or the ampli¬ 
tude of “roughness”) er of a membrane is of the order of 
the inter-membrane separation d, viz., a 2 = fid 2 , where 
ft is some numerical constant. One may interpret the 
mean square fluctuation as the inverse curvature of a 
harmonic potential well. By integrating out the bending 
fluctuations and applying a 2 = jld 2 , one then obtains a 
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steric potential of the form Vjj = cufk^T) 2 /nd, 2 , which 
is known as the Helfrich fluctuation potential HE], with 
k the curvature stiffness of the membrane. From this po¬ 
tential one derives the fluctuation osmotic pressure, one 
that is entirely generated by constrained thermal fluctu¬ 
ations of the membranes. It has been argued (see, e.g., 
EH3E3]) in the past that eg is a universal number 1 141 . 
On the other hand, one may reasonably expect that the 
rms fluctuation a, being the effect of thermally activated 
undulations of the membrane, vanishes as T —> 0 and/or 
k —i oo , and accordingly fi should depend on T/k. More¬ 
over, as d increases, the behavior of a should cross over to 
that of a freely undulating membrane, i.e., scale as \/S 
(where S is the transverse projected area of the mem¬ 
brane; cf. Ref. m, and thus a should also be a (possibly 
nonlinear) function of d/V~S. 

In the foregoing paragraph we have described the prob¬ 
lem of the steric potential generated purely by the ther¬ 
mal undulations of a pair of membranes and the hard-wall 
constraint. In more realistic systems (including those of 
biological interest), the steric potential is modified by 
the presence of other interactions, both short- and long- 
range, and these may include van der Waals, electrostatic 
(if the membranes are charged), charge regulation and/or 
hydration forces. Problems of relevant concern include: 
(i) what is the behavior or form of the effective inter¬ 
action potential? (ii) How does the equilibrium inter¬ 
membrane separation change as a function of external 
osmotic pressure? (iii) In the case where the interaction 
is attractive at long range but repulsive at short range, 
is there an unbinding transition at zero external osmotic 
pressure? If so, what is the order of the unbinding tran¬ 
sition and the threshold value of the (attractive) inter- 
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action strength for the membranes to unbind? Besides 
the presence of additional interactions in real systems, 
the membranes themselves can be multi-component, e.g., 
consist of a mixture of lipids and lateral inclusions that 
may undergo phase separation EE and the solution may 
also contain polymers (such as proteins) which then in¬ 
teract with the membranes and anchor or adsorb onto 
the membrane surfaces Ell- 

Theoretical treatments of such problems have tradi¬ 
tionally fallen into four classes: mean-field or Flory- 
Huggins-type approaches (see, e.g., E3 EMU) , vari¬ 
ational Gaussian approximation (VGA) schemes (see, 
e.g. , [3H7J E]), FRG-based methods (see, e.g., [T2], [501 - 
and Monte Carlo simulations of different flavours 
(see, e.g., [2M33] ). For the simpler case of a pair of 
single-component membranes interacting via a direct po¬ 
tential that consists of a long-range attractive van der 
Waals tail and a short-range repulsive hydration poten¬ 
tial, mean-field theory (MFT) based on the simple addi¬ 
tion of the direct potential with the fluctuation-induced 
steric potential predicts a first order unbinding transi¬ 
tion at a critical strength of the Hamaker coefficient EE 
whereas a Flory-Huggins-type theory, modeled after the 
van der Waals theory of liquid-gas condensation, predicts 
a continuous unbinding transition at a critical strength 
of the Hamaker coefficient W = W c and scaling behav¬ 
ior d ~ | W — W c | _1 [TS]. Such an unbinding transition 
in three spatial dimensions has also been predicted by 
FRG-based calculations [20) 21]. On the other hand, ex¬ 
perimental reports that the unbinding transition in lipid 
bilayer systems is of first order (see e.g. Refs. [54H3T] f 
have been more than those that report the unbinding 
transition to be of second order (see e.g. Ref. [38]). 

Using a novel approach based on a self-consistent treat¬ 
ment of the hard-wall constraint, we revisit the prob¬ 
lems described in the foregoing paragraphs. The require¬ 
ment of self-consistency is realized by implementing the 
hard-wall constraint at the level of the partition func¬ 
tion via the use of Heaviside function. By making use 
of a representation first proposed by Panyukov and Ra¬ 
bin [35] in a rather different context, the steric constraint 
can be transformed into terms of an effective interac¬ 
tion Hamiltonian. By considering small undulations of 
the membrane, we derive closed-form self-consistent ex¬ 
pressions for the mean square fluctuation er 2 and the 
equilibrium inter-membrane separation distance, which 
precisely describe how the strength of fluctuations is re¬ 
lated to both temperature T and the curvature stiffness 
k of the fluctuating membrane |T], as well as the aver¬ 
age inter-membrane separation distance d. This effective 
Hamiltonian also enables one to see the precise mecha¬ 
nism in which the fluctuation osmotic pressure emerges 
from fluctuation and steric forces, and shows a crossover 
of scaling behavior from a 1 ~ \Jk-QT/ nd 2 , when a pair of 
membranes are close to one another, to er 2 ~ {k-sT/rijS 
(where S denotes the transverse projected area of each 
membrane), when the membranes are far apart. 

Our self-consistency requirement is reinforced for the 


case of real systems (i.e., where the direct potential 
now includes non-Gaussian terms stemming from elec¬ 
trostatic, van der Waals, charge regulation, and/or hy¬ 
dration interactions) by our use of the VGA. The VGA 
captures the fact that the magnitude of the fluctuation- 
induced steric potential is influenced by the direct po¬ 
tential, and vice versa U El SO]- Mathematically this 
translates into the problem of solving a coupled pair of 
(typically nonlinear) equations for the rms fluctuation 
amplitude and the average separation (or more generally 
speaking the fluctuation correlator of an observable and 
the average value of that observable). It is known that 
the VGA becomes exact as the co-dimension of the man¬ 
ifold (which is equal to the number of components of the 
vector field representing the transverse fluctuations of the 
manifold) tends to infinity |41j . 

The VGA can be regarded as being complementary 
to FRG techniques. The FRG is essentially an asymp¬ 
totic method which has been used to study problems 
relevant to critical behavior and the associated scaling 
exponents (i.e., the “shape” of the effective interaction), 
addressing the limiting behavior of the interaction po¬ 
tential in the regions of small and large inter-membrane 
separations. As such it is less useful as a tool for making 
predictions about the system’s behavior at intermediate 
separations, which are the typical length scales of bio¬ 
logically relevant systems. At such length scales there is 
usually a number of competing interactions of compara¬ 
ble strength, and the use of asymptotic methods is not of 
much practical value [52]. One of the primary interests 
in membrane biophysics is also to elucidate how the ef¬ 
fective mesoscale inter-membrane interaction, accessible 
via detailed osmotic stress experiments on multilamellar 
lipid systems [401 I43H5B] , arises from details of more mi¬ 
croscopic forces, an aspect that is not captured by the 
FRG-based methods at all. These aspects can be ade¬ 
quately addressed by a VGA-based approach. Moreover, 
the VGA-based approach is comparatively simple to im¬ 
plement and the physical mechanism shines through the 
formalism more transparently. 

Among our main motivations here is also to analyze 
the behavior of the electrostatic correlation interaction 
in the strong coupling limit m that has been already 
observed between stiff charged silica surfaces @2 and is 
now being investigated also in the case of soft, fluctuat¬ 
ing membranes [45j . It seems to us important to be able 
to understand the difference between the experimentally 
measured ion correlation effect between stiff and soft in¬ 
terfaces and to isolate in what way the membrane fluc¬ 
tuations change the expected separation behavior of the 
strong coupling interactions. Our analysis is the first in 
this direction and should be helpful to understand the 
experimental data when available. 

Our Paper is divided into the following sections. Sec¬ 
tion II introduces the problem of determining the os¬ 
motic pressure of a freely undulating membrane near a 
hard wall, and describes how the hard-wall constraint 
can be transformed into effective potential terms via the 
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Panyukov-Rabin representation. The system considered 
is a pair of membranes subject only to hard-wall repulsion 
and undergoing thermal bending fluctuations. The first 
main result of this section is Eq. (2.23) for the rms fluctu¬ 
ation amplitude of a membrane near a hard wall, which 
interpolates between the scaling behavior of a confined 
membrane at short separations and that of a freely undu¬ 
lating membrane at large separations. The second main 
result is Eq. (2.311, which expresses the steric potential as 
the sum of contributions from thermal bending and zero 
mode fluctuations. In Sec. Ill, we introduce a variational 
method due to Feynman and Kleinert j50| , which we ap¬ 
ply to study the behavior of membranes in more realistic 
systems, i.e., ones that are characterized by effective in¬ 
teraction potentials of non-Gaussian form (specifically, 
an attractive square well potential, the Morse potential, 
the Moreira-Netz potential, and a combination of van der 
Waals and hydration energies). In Sec IV we present our 
discussion and conclusions. 


II. MEMBRANE NEAR A HARD WALL 
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FIG. 1: Comparison of plots of 0 a(*) generated using the ex¬ 
act function [Eq. (2.91; behavior is described by the blue-grey 
circles, yellow squares, and green diamonds for A = 0.01a, 
0.1a and 0.3a respectively, where a is a microscopic length- 
scale] and the saddle point approximation [Eq. ( |2.13 1; behav¬ 
ior is displayed by the blue, red dashed, and green dot-dashed 
lines for A = 0.01a, 0.1a and 0.3 a respectively]. The agree¬ 
ment is quite good, with the accuracy improving for smaller 
values of A. 


The fluctuation osmotic pressure of a pair of mutually 
impenetrable membranes of bending stiffnesses «q and 
pL] is equivalent to that of a membrane of stiffness k = 
kik 2/(^1 + K 2 ) near a hard (i.e., impenetrable) wall. We 
shall thus consider a single membrane fluctuating near a 
hard wall as our prototypical system. To maintain the 
membrane at constant average distance from the wall, an 
external osmotic pressure P that is equal and opposite to 
its fluctuation osmotic pressure has to be applied. We fix 
the position of the wall at the origin of the z-axis, and a 
point on the membrane shall have transverse coordinates 
xi = (x,y) and occupy a position f(xj_) on the z-axis. 
The space between the wall and the membrane could be 
empty or filled with electrolyte solution. The membrane’s 
total energy is given by 

H = |dV(^(V i£(xy)) 2 + VWxy))), (2.1) 

where the interaction potential is given by 

V{l{x±)) =w(£(x ± ))+Pi(x ± ). (2.2) 


integral over the transverse projected area of the mem¬ 
brane) , and the resulting integration over the fluctuating 
fields is trivially Gaussian. 

There are three terms in Eqs. (2.1) and \2.2\ . The 
first is the free energy associated with bending undulation 
inodes of the membrane. The second term w represents 
the direct potential between the wall and the membrane, 
arising from forces of electrostatic and non-electrostatic 
(such as hydration) origin. In general the direct poten¬ 
tial has a non-trivial form and is non-Gaussian in fluc¬ 
tuations. The third term represents the effect of apply¬ 
ing an external osmotic pressure P on the fluctuating 
membrane. Experimentally, for membranes in solution 
the external osmotic pressure can be controlled by using 
the classic method of Rand and Parsegian [33], which 
is to vary the concentration of polymer (such as dex- 
tran |44j . polyvinylpyrrolidone |4Sj or polyethylene gly¬ 
col [46]) added to extract water from between the mem¬ 
branes. Theoretically, the osmotic pressure term can also 
be interpreted as a source field that one uses to differenti¬ 
ate the logarithm of Z in order to obtain the equilibrium 
separation distance (£q). 


Denoting the (geometric) mean separation by £n, defined . TT . „ , . . n ,. . . 

1 r 79 " A. Hard-wall constraint: lunctional representation 

by J d z x± t(x_|J = t 0 , we can write 


*(x±) = £ 0 + 6£{x ± ) = 4 + J^ e^Sig, (2.3) 


where the wave-vector integral excludes the Q = 0 mode, 
and 51q are the F ourier modes of the deviation 6£(x _l). 
Very soon [cf. Eq. (2.21)] we shall see the beauty of such 
a decomposition that lies at the heart of the Feynman- 
Kleinert variational theory: the terms that are linear in 
fluctuation vanish from the Hamiltonian (owing to the 


The partition function of a fluid membrane near a hard 
wall is then given by 


{x-l} 


/ oo _ 

<^(x_l) 0(f(xj_)) 

-OO 


(2.4) 


where j3 = l/k^T and 0(K(xjJ) is the Heaviside function 
which is equal to unity if £(xj_) > 0 and zero if ^(xj_) < 0. 
Taking up an idea in Ref. [33, we express the Heaviside 
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function in terms of the Dirac delta function: 


B. Steric potential: saddle-point approximation 


/»oo 

0(£(x x ))=/ dAS(A-i(x ± )). (2.5) 

Jo 


Next, we use a Gaussian representation for the delta func¬ 
tion, and write the Heaviside function as 


e(«(x±)) = 


lim 0 a (£(xjJ) 

A—>0 


= lim / dA 

A—^0 ,/n 


(A-l(xjJ ) 2 


a/27tA 2 


( 2 . 6 ) 


The limit A —> 0 is a mathematical idealization for an 
“infinitely hard” wall. However, for real membranes, the 
boundary region is not infinitely sharp, but has a certain 
extent which is set by the lengthscale of the lipid head- 
group size a. The parameter A reflects the size of this 
boundary region. We can thus write 


A = 


(2.7) 


where c depends on the microscopic details of the chem¬ 
ical make-up of the membrane. 

By making a change of variables A = f u 2 and pro¬ 
moting a prefactor u into the exponent by means of a 
logarithm, Eq. (2.6 1 can be expressed in the form 


/ OO 

dn 

-OO 


1 

a/27tA 2 


of\(. u A x ±)) 


where 




2A 2 


( 2 . 8 ) 


(2.9) 


We have obtained an exact functional integral represen¬ 
tation of the steric constraint for each point of the mem¬ 
brane. In contrast to A, the variable u now runs over 


the unbounded interval (— 00 , 00 ). Using Eq. (2.8), we 


can express the partition function in Eq. (2.4) as 

O 

d£(x.±) 


Z = lim IT 
A->0 J -- L 
bi}' 

_ r d 2 x ± (A(x ± )-!(x ± )) 
P J T2 o.\2 


dA(x x ) 


0 a/27tA 2 

~PH 


= iTo n / di ^ 
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r°° du{x_ l) 
Loo a/ 27 tA 2 


J- 


f A(«(xj_)b(xj.)) e ~0H 


( 2 . 10 ) 


In the functional integral expression for Z , we have in¬ 
troduced a microscopic cut-off length scale a, in order to 
convert the sum over a two-dimensional lattice of points 
into an integral. The length scale a is set by the size of 
the membrane’s molecular constituents, e.g., the diame¬ 
ter of the lipid headgroup. 


Let us now make the saddle point approximation to 
e^ x . It is given by 




27 T 


l/"K P T)l 


0 A(“s pb) 


( 2 . 11 ) 


where rt sp solves the saddle point equation 

/aKpT) = 0 


sp 


(Xj.) = Tx x ) + vL(xl) 2 + 2A 2 (2.12) 


We find that the saddle point approximation to f\(u) is 
given by , where 


/a = 


2 2 | A 2 T + a/ 2 A 2 + £ 2 ) j 

+ * \n{l + x/2A 2 + £ 2 ) - ^(* ~ V2X 2 +i 2 ) 2 , 

(2.13) 


where we have written £ = £(x x ) to lighten our notation. 
In Fig. [l] we compare the saddle point appro xima tion of 
Q\ (viz., e-L) with the exact function [Eq. (2.9)]. The 
agreement is excellent, especially for smaller values of A. 
We can thus approximate Z in Eq. ( 2.10| ) by its saddle- 
point value, viz., 


Z* = lim TT 
) 

bil' 


A->0 


d£(x _l) 


J - 


Ii/lbbi)) e -P h 


(2.14) 


On going from Eq. (2.4) to Eq. (2.14), we have effectively 


moved from a partition function involving £{x±) subject 
to a hard-wall constraint to a partition function where 
£(x_l) is unconstrained , but the steric condition is en¬ 
forced energetically through terms in an effective Hamil¬ 
tonian. In our calculational steps we have not made any 
assumptions and made use of only one approximation, 
viz., the saddle-point approximation of Eq. (2.11). 


As A is small, we can expand the above to leading order 
in A: 


/a = 


P + eVF 1 vT 2 i, 

-4 + TT + 2 ln 


4 A 2 


e + Vp 


2 i 


< f + 5 ^))L + lllA + 0(A 3) 


+ Vf 2 ) 


(2.15) 


where we have neglected a constant term |ln(27r). We 
see that terms for which £ is negative diverge to negative 
infinity, and so e-G goes to zero, which reflects the zero 
probability of finding the membrane inside the wall. 

We now focus on the case of membranes undergoing 
thermally excited undulations with small amplitude. Let 
us write £(x x ) = £0 + 8£{x j_), where £q is the geometric 
(not thermal) mean defined by £q = {1/S) f d 2 x±£(x±), 
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and <5£(x_l) are small deviations around 4- In Fourier 
space, 4 is the zero wave-vector mode and 

W(x_l) = J2' eiQ ' X±S£ Q (2.1.6) 

Q 


where Q is a two-dimensional wave-vector conjugate to 
xj_, and the prime denotes the exclusion of the zero wave- 
vector mode from the wave-vector sum. Both 4 and 
{5£q} [or W(xj_)] are independent thermally fluctuating 
variables that take values from the interval (— 00 , 00 ) [cf. 
Eq. (2.10l]. By performing an expansion to quadratic 
order in 5£{x j_), Eq. (|2.15[) becomes 


(4 + 5y^g)A 2 (lo + Sy/QA 2 §£ 
8£oVP 0 (Zo + V?o) 4(£ 2 )3/2(4 +y^) 
(24 + y/lf)(4 + 5y^f)A 2 ^ 2 ^ 

+ V^o) 2 


For negative values of 4 the leading-order term of the 
above expression diverges to negative infinity, which 
means that the geometric average position of the mem¬ 
brane cannot move inside the wall. Under this approxi¬ 
mation, we can effectively restrict the range of values that 
£0 takes in the partition function to the positive interval 


(0, 00 ), and re-express Eq. (2.17) as 


/; 


3A 2 3A 2 


We note that the term linear in <W(x_i_) will vanish when 
we integrate over the transverse projected area. The sad¬ 
dle point approximation to Z can thus be written as 


r uvj r 

Zxs d£ 0 U d6£{x±)t 


,-0fd 2 x x (*(V 2 x i(x x )) 2 +V) 


{x-l} 


where V is now an effective potential defined by 


(2.19) 


V = w(£(x±)) + P£(x ± ) (2.20) 

and 

... .. .. 9k B Tc 2 2 3 k B Tc 2 

w(£(x ± )) = w(£{x_ l)) + ffl(xi)) 2 + 8£2 . 

(2-21) 

In the above, we have removed a prefactor (27rA 2 ) _s / 2a . 
This can be done, as the prefactor simply corresponds to 
a shift of the effective free energy by an (infinite) con¬ 
stant. 

The prefactor of the term quadratic in S£ can be in¬ 
terpreted as an effective compression modulus gener¬ 
ated by the steric force acting between the membrane 
and the wall. The interaction potential w thus consists 
of two contributions: a soft potential w that describes 
the longer-range (and coarse-grained) interactions of the 
membranes, and a hard potential which arises from the 
steric constraint. Such a decomposition into soft and 


hard contributions is possible because we are working 
at an already coarse-grained/mesocopic level; otherwise 
all contributions (including steric effects) are reducible 
to electrostatic ones. By rewriting the steric constraint 
in terms of delta functions [Eqs. (2.5) and 2.8)] and ap¬ 
plying a saddle point approximation on u, we were able 
to transform the partition function for a sterically con¬ 
strained membrane into one for an unconstrained mem¬ 
brane, with the effects of the steric constraint accounted 
for by terms in an effective interaction potential w [cf. 
Eq. (2.21)]. This also completes the program set out in 
Ref. [S] but not carried out to its full implementation. 

Two observations can immediately be made: (i) the 
energy diverges as the mean separation 4 tends to zero, 
which reflects the fact that the membrane is unable to 
penetrate the wall, and (ii) the prefactor of the fluctu¬ 
ation term also diverges as £0 —> 0, which reflects the 
fact that the fluctuations of the membrane must also be 
suppressed as it reaches the wall. 


C. Application: steric and fluctuation forces 


Let us consider a pair of membranes that interact via 
only steric and fluctuation forces. As we have already 
noted, the two-membrane system can be recast as a mem¬ 
brane interacting with a hard wall. This problem admits 
of a closed-form solution that describes how a 2 depends 
on 4- We set the soft potential to zero: w = 0. The hard- 
wall and fluctuation-induced repulsion from the wall has 
to be balanced by an external osmotic pressure if the 
membrane is to remain at a finite distance from the wall. 
For the potential including only the sterically generated 
interactions, Eq. (2.21) assumes the form 


w(£(x±)) = 


9k B Tc 2 

~W~ 


0«(xi)) 2 


3 k B Tc 2 


( 2 . 22 ) 


The functional integration over 5£ in Z is thus Gaussian 


and can readily be performed. Using Eq. (2.19) and the 


definition cr 2 = ((<W(x_i_)) 2 } st (where (...)« denotes aver¬ 
aging over 5£ using Z for a given 4) yields for the mean 
square fluctuation amplitude [cf. also Eq. (3.9)] 


<x 2 = 




4 

(o_ 

12c 


7 r tan \3cS\ k B T 


(2.23) 

Comparing with the Ansatz a 2 = fid 2 of Ref. [2] (where 
d = 4)> we see that fi has the following structure: 


V = 


12c 


—(l — — tan 
k \ 7r 


-1 


2£ 2 


3cS A/ k B T 


(2.24) 


In accordance with the expectations described in the 
Introduction, our calculation has revealed the structure 
behind the prefactor fi, showing it to be a nonlinear func¬ 
tion of k B T /k and £q/\[S. To understand the properties 
of fi, we look at two limiting cases: the regimes of small 
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and large inter-membrane separations. Defining the di¬ 
mensionless variables £q = £o/d* and d = cr/d*, where 
d* = v / 3c(fc B T) 1 / 4 v^/(4K) 1/ ' 4 , the above equation can 
be put in dimensionless form: 

( “ 5) 


To order S 3 £ 0 s , f s is given by 

. p/ , 3k ^ T ° 2 , 9 (k B T) 2 c 2 S 27 (k B T) 3 c 4 S 3 

s 0 8£q 32ttk£q 1287 r «: 2 t'Q 

(2.30) 

In the near-field regime, £ 0 -C \/3c(k B T/4:K) 1 ^ 4 \ /r S. To 
order S’ -1 , f s is given by 


At separations small compared with the linear dimen¬ 
sion of the membrane, viz., £ 0 < y/3c(k B T/ 4k) 1 / 4 -\/<S', we 
obtain 


a 


2 


1 

12 c 



(2.26) 


This formula has the same scaling dependence on sepa¬ 
ration as the one originally postulated by Helfrich, but 
here we have derived the scaling dependence rather than 
postulated it. From our result we see that the fluctu¬ 
ations become small at low temperature and/or large 
membrane curvature modulus. On the other hand, for 
two membranes that are more widely separated apart 
than the linear dimension of either membrane, viz., £ 0 > 
\i3c{k B T/An) 1 P\/S 1 there is a crossover to the behavior 
of a single, free membrane: 


2 _ k B TS 

47 TK 


(2.27) 


At large separations, each membrane would behave as 
if there is no hard wall potential present, and the mean 
square fluctuation of each membrane is then set by its 
total area. 

The free energy per unit area f s for a given average 
separation £q and external osmotic pressure P is given 
by fs = —T In Z. The functional integration over the 
fluctuation modes 5£q in Fourier space is Gaussian and 
yields 


fs — P£ o + 


3 k B Tc 2 


k B T £ d 2 Q r 9k B Tc 2 \ 

2 7(277)2 +4 k£ 4 Q 4 }' 


(2.28) 

In the above calculation, we have subtracted off a con¬ 
stant background contribution to f s coming from £q —> 
oo eh- Taking the upper bound of the momentum inte¬ 
gral to be oo and the lower bound to be 1/v S, we obtain 
a closed-form expression for f s : 


fs — P£ 0 + 


3 k B Tc 2 


k B T 
87 tS 


In 


9k B Tc 2 S 2 \ 
4k£ a J 


, k B T 1 9k B Tc 2 f ^ V2(4 K y/ 4 £ 0 ) 

+ 4tt y 4 k£ 4 l an V {9k B Tc 2 ) l / 4 \/~S ) 


- tan 


-1 


x v^(4k)V% \\ 
(9k B Tc 2 ) l / 4 VSj 


(2.29) 


We can consider the far- and near-field behavior of 
f s . The far-held regime is described by £q 


3 k B Tc 

k B T 
4t tS 


kBT + -k^ + P£o 




1 + ln 


/ 3 C\fki 


TS 


V 2,/k^o 


(2.31) 


Equations (2.301 and (2.311 indicate that the steric po¬ 


tential has two contributions: one that comes from ther¬ 
mal fluctuations of the zer o mod e (i.e., £q), r epres ented 
by the second term in Eq. (2.301 and in Eq. (2.311, and 


another that is induced by thermal bending fluctuations 
of the membrane (i.e., <W(x_l)), represented by the third 
and fourth terms in Eq. (2.30) and the first and last terms 


0 to £ 0 4 , and 


in Eq. (2.31). Note that as £0 is varied from small to large 
values, the bending fluctuation contribution to the steric 
potential changes its scaling form from 
is thus much weaker at large separations. 

In the near-held regime, to zeroth order in S~ 4 , we can 
compare our result above with the picture of Ref. [2], 
where the free energy cost of thermal bending fluctu¬ 
ations of steric membranes is described via the term 
Vh = Cf{(k B T ) 2 / k£q\ in keeping with the literature we 
call Vh the Helfrich interaction. In the Helfrich inter¬ 
action term, ca is regarded as a universal number (see, 
e.g., Refs. (2j [13]). One may regard the first term of 


Eq. ( 2.31[ ) as being the analogue of the Helfrich interac¬ 
tion 49]. However, one should note an important distinc¬ 
tion: for k, —y 00 , Vh goes to zero, and the Helfrich term 
thus does not address the case of osmotic pressure gen¬ 
erated by longitudinally fluctuating hat membranes (i.e., 
zero mode huctuations) and therefore describes only the 
bending huctuation “decoration” about a hat membrane, 
but not the longitudinal huctuation of the hat membrane 
itself. 

As the intervening region between the membranes is 
free to exchange volume with the surrounding reservoir 
whilst being subject to a constant external osmotic pres¬ 
sure P , the system behavior is best described using a 
constant osmotic pressure ensemble, where the equilib¬ 
rium separation (£q) is determined according to 


1 f c 

w s zL 


d£n £1 


0 e 


-psf. 


TdlnZ 
S f)P 


(2.32) 


The equilibrium separation can be estimated in the 
saddle-point approximation by looking for the mini¬ 
mum of f s with respect to variations in £q. Defining 
the dimensionless variables £q = £o/d*, p = (3P(d*) 3 , 
S = S/(d*) 2 and a = 3c\/k B T / 87 Tywhere d* = 
\/3c(k B T) 1 ^ 4 \/S/(An ) 1 / 4 , we can recast Eq. (2.29) in di- 
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FIG. 2: Steric interaction: behavior of the dimensionless sep¬ 
aration £* as a function of dimensionless osmotic pressure p 
and (inset) behavior of dimensionless rms fluctuation ampli¬ 
tude a as a function of £*, where £* « {£o)/d* is the equilib¬ 
rium separation in the saddle-point approximation, a = cr/d *, 
p = l3P(d*) 3 , and d* = V3^(A;bT) 1/4 V^/(4k) 1/4 . 


mensionless form: 

Is = (3{d*) 2 f s ( 2 . 33 ) 

= o + ( 3 c 2 / 8)£ 0 2 — — In (l + £ 0 4 ) 

+a£g 2 (tan _1 (l — V2 £ 0 ) + tan _1 (l + V^f'o)) 

The stationarity condition df s /d£ 0 = 0 then leads to 

_ ^ + 2a( tan _1 (l - V2£*) + tan" 1 ^ + y/2£*)) 

P = ---~-, 

(P ) 3 

( 2 . 34 ) 

where £* is the saddle-point approximation to the equi¬ 
librium separation (£q). In the near-held limit (£* <C 1), 
we obtain to leading order 


The behaviors of £* and a are plotted in Fig. [ 2 ] 
The rms fluctuation amplitude vanishes at zero inter¬ 
membrane separation (reflecting the hard-wall con¬ 
straint), scales linearly with the inter-membrane sepa¬ 
ration for small values of the separation, and saturates 
at a constant value set by the cross-sectional area of the 
membrane for large values of the separation. 


III. FEYNMAN-KLEINERT VARIATIONAL 
APPROXIMATION 


In the previous section we derived an effective Hamilto¬ 


nian [Eq. (2.21)] which accounts for the steric constraint 
via the hard potential, and applied it to study a mem¬ 
brane interacting with a hard wall via only Helfrich and 
steric forces. We were able to obtain closed-form ex¬ 
pressions for the rms fluctuation amplitude and free en¬ 
ergy per unit area of the membrane [cf. Eqs. ( 2.23[ ) and 
(2.29)] because for such interactions the form of the ef¬ 
fective Hamiltonian is Gaussian. In realistic systems the 
soft potential rarely has a Gaussian form, and further ap¬ 
proximations will have to be made on the partition func¬ 
tion. One such approximation is the Feynman-Kleinert 
(FK) variational approximation |5(’)j , developed originally 
for the quantum-mechanical partition function of an an- 
harmonic oscillator, but can equally well be applied to 
the classical partition function of thermally fluctuating 
membranes 

The application of the variational approximation [52] 
begins with a trial partition function Z\ that is Gaussian 
in fluctuations 8£, viz., 


2i = 


/»OO f OO / \ 

d £ 0 1 ] / d «( x J .) e -^ d2 ^(^ v ^( x -» 2+y ‘-0 

J 0 r__ 1 J — OO 

(3.1) 


I*-!-} 


where 


/?* \ — 3 


(2.35) 


Vtrial = \b («( Xj J ) 2 + w 0 (£ 0 ) +P£{xl±). (3.2) 


Restoring dimensions, we have 

W .(W)V3(^(I^ +c )) 1,3 p-V3. 

(2.36) 

Thus, the equilibrium separation diverges as P -1 / 3 as 
P —> 0, and vanishes as T —> 0 (the pair of membranes 
simply collapse onto each other at zero temperature) as 
one would expect. 

In the far-field regime (£* 1), we have to the order 

of (P )" 5 

+ WW* ( 2 .37) 

4(P ) 3 47t(P) 5 


The variational principle relies on Jensens’ inequality, 
viz., 

Z>e - (3.3) 

where the notation (.. .)i denotes Boltzmann averaging 
with respect to the statistics of both £q and {<5£(xj_)} 
specified by Z\. The best estimate is obtained by looking 
for the maximum upper bound on the right hand side. 
There are two unknowns: (i) an unknown parameter P, 
which is related to the mean square fluctuation of the 
membrane, and (ii) an unknown function wq. The form 
of Wq is fixed in terms of B by requiring that 

(V - Rtrialh = 0 . (3.4) 


The equilibrium separation increases monotonically with The unknown parameter B is then determined by opti- 
decreasing osmotic pressure. mizing Z\. 























and the associated average (... ) a , viz., 


To see how the FK variational approximation works 
in the context of fluctuating membranes, let us define a 
restricted trial partition function Z a , viz., 





dS£(x ± )e~^ d2 ^ < f (Vl<(*-L)) a +*W) ; 

(3.5) 




n {xi} fZo dS£{x_ l) (...) e~PS d2x ^ f ( V i^( x r)) 2 +Vtriai) 

n {x } /^ 0 dW(x_ L )e“^/ d 2 a: - L (4( v i^ x - L )) 2+Vtrial ) 


(3.6) 


By using Eqs. (2.201, (3.2 1 and (3.4), we obtain 


/ OO 

dM(x x ) (F - F tria i) = o 

-oo 

IX-i-J 


wo(£o) = w a 2 ( 4 ) - — a 2 , 


(3.7) 


where the variational estimates of the mean square fluctuation a 2 and the interaction energy w a 2 are defined by 


a 2 = ((<W(x x )) 2 ) a ; 
«V 2 (4) = (w(£))a- 


(3.8a) 

(3.8b) 


We compute the mean square fluctuation: 


a 2 = 


[<?x± ? n {Q > 0} f dSe%dS£$((S£$) 2 + (^ m ) 2 )e 

J s X± -p S- 1 ^ OQ 4 + 

n {Q >o}/^Q^ n e 


-/SS- 1 ^ ( K Q 4 +B)((«g) 2 +(«“) 2 ) 


Q>0 


Q>0 


+B)((« r Q G 2 +(«g 1 ) 2 ) 


^ (l-*ten-'(U£ 


8 VkB V ^ 


S \ B 


(3.9) 


In the above calculation, we have made use of the con¬ 
tinuum representation of the wave-vector sum, S ' -1 )T) —» 

{Q 1 

/ ( 2 ^ , and restricted the wave-vector sum to those wave- 
vectors that are positive, the reason being that <5£(xj_) is 
real and thus the components S£q (and S£ 1 q) are not 
independent: 5£q = S£ T fiq and 5£'q = — Sifiq . In the 


second line of Eq. (3.9), we have included a correction 


term due to the finite size of the membrane. In the limit 
of an infinitely large membrane (which effectively means 
that the square root of the membrane’s projected cross- 
sectional area is much greater than the inter-membrane 
separation), the above result simplifies to 


a 2 = 


k B T 
8 \/kB 


and 


B = 


(k B T) 2 

64ft(X 4 


(3.10) 


(3.11) 


which agrees with Eq. (10) in Ref. [ 6 ]. Let us compute 
w a 2 (£q) by taking explicitly into account the Gaussian 
variational Ansatz , leading to 

/ °° dk 

—w(fc)(e lW(x - L) ) CT 


— f d£w(£) 

-OO 

dl //1X _ (i-io ) 2 


/ 27 T 

r —00 J —OO 


/-ooV27TCr 2 


w(£) e 


(3.12) 


In the second line, we have performed a second order cu- 
mulant expansion for the Gaussian variational Ansatz , 
followed by an inverse Fourier transform of w, and £ is a 
dummy variable [not to be confused with the thermally 
fluctuating field £(x x ) of the first line] that runs from 
—00 to 00 . Next, we have integrated over k. The fi¬ 
nal result thus depends only on £q and a. Knowing the 
form of the interaction potential w(£), we can plug it into 
Eq. (3.12) to obtain its variational estimate. 
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In what follows, it is useful to define the variational 
free energy for a given £ q: 






(3.13) 


To determine the variational free energy, we first express 
Eqs. (3.1) and (3.2) in Fourier space: 


Z-, = Jd £o H J(IMqJ dd£^ l e~ fis ^ 0+Wo{eo)) 

(3.14) 


(Q>oR 

- Ps - 1 E (« Q 4 + b )|« q | 2 


xe 


Q> o 


Integrating over real and imaginary modes of S£q, we 
obtain 


On the other hand, / var is not applicable to the strong 
fluctuation regime, defined to be one in which the fluc¬ 
tuation potential is stronger than the attractive tail of 
the direct potential, because in such a regime the fluctu¬ 
ations are strongly nonlinear and the Gaussian approxi¬ 
mation assumed by the VGA is no longer reliable. The 
VGA-based / var thus applies to a regime intermediate be¬ 
tween the weak and strong fluctuation regimes, which as 
already stated, is exactly the regime we want to address. 

To obtain a relation between er 2 and £q, we vary / var 
with respect to a 2 ; this yields 


dw a 2 (k B T) 2 9k B Tc 2 

da 2 128kct 4 + 8Kq 

This can be rewritten in the form 


(3.20) 


z .= n 

{Q>o} J 
/ nk B T 


f . -S - 1 E (*Q 4 +b) |« q | 2 

d6£% I ddi'^e «>° 


n 

{Q>o} 


\ kQ* + B 




= e Q>0 


' »!' ,,ir 'A > yry) = e ~fy/£ (3. 


15) 


The regularisation subtracts off the constant (and diver¬ 
gent) contribution of membranes that are infinitely far 
apart 155 ■ Equation (|3.14[) thus becomes 


Z 1= \dt oe -^(« 0 +»°(W+¥v^). 


(3.16) 


Subst itutin g for wq(£q ) its value from Eq. (3.7) and using 
Eqs. (2.21) and ( |3.8b ), we obtain 


Z 1 = dioe 


(3.17) 


or equivalently, from Eqs. (3.10) and (3.11) 


Zi = / dJU 


oe 


~/3Slw 2+ (fe B'T) + 9fc B Tc <r 3k B Tc p£ 

n r V 128 KCJ 2 8 £% 8 & 


. (3J8) 

Comparing with Eq. (3.13), we see that the variational 


free energy per unit area is given by 

f _ M 2 , 9 k B Tc 2 a 2 , 3k B Tc 2 „ ine /o 
/var - 1 r,o.._9 H o~/)4 -1-^2- \-W a 2+Pl Q. (3.19) 


128k<t 2 


8£ 2 


In the language of Ref. JB], the first three terms can be 
regarded as comprising a repulsive fluctuation potential. 
The fourth term is a fluctuation-renormalized version of 
the direct potential, and thus / var cannot be simply re¬ 
garded as a sum of fluctuation and direct potentials. The 
additive sum would be valid in the so-called weak fluctu¬ 
ation regime, defined by Ref. [26] to be one in which the 
attractive tail of the direct potential is stronger than the 
repulsive fluctuation potential. 


dw a 2 1 9k B Tc 2 
da 2 = 2 B 8£ 4 ’ 


(3.21) 


which is similar to Eq. (17) of Ref. [6], with the extra 
contribution coming from the steric potential. 

The equilibrium value of £q for a given external osmotic 
pressure P is defined by 


(4)i 


l 

MP) 



(3.22) 


The equilibrium separation (£o)i is thus a function of P. 
We invert the above relation to determine how P depends 
on (£q)i, which gives the equation of state. For this one 
generally has to resort to numerical means. On the other 
hand, one can approximate (£o)i by its saddle-point value 
£q, obtained by minimizing / var over £q\ 


'//var 

d£o t 


(3.23) 


This relation yields £q as a function of the external os¬ 
motic pressure P. To summarize: in the Feynman- 
Kleinert approach, one decomposes the thermally fluc¬ 
tuating field into two types of contributions, a zero¬ 
mode contribution £o and a finite wave-vector contribu¬ 
tion {^q}q/oi approximates a 2 by means of a Gaussian 
kernel but retains the full nonlinear dependence of the 
free energy on £q (present in w a 2 ), which is then mini¬ 
mized over £q. This leads to an improved accuracy over 
that of the more conventional form of VGA [S2], which 
approximates both the thermal equilibrium average of an 
observable and its mean square fluctuation by means of 
a Gaussian kernel. 


A. Applications 

In what follows, we apply our formalism and the vari¬ 
ational approach developed in previous sections to study 
the physical behavior of four different model membrane 
systems in the regime where £q -C y/S. 
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First, we analyze a system consisting of equally 
charged impenetrable membranes, with an intervening 
solution of multivalent counterions, assuming that the 
surface charge density is sufficiently large so that the 
system is in the strong coupling (SC) regime of elec¬ 
trostatics [53H56] . In this regime, for the case of two 
planar, infinitely rigid membranes, the osmotic pressure 
decomposes into a simple sum of two contributions: one 
describing the entropy reduction of each counterion and 
another describing the electrostatic interaction between 
the counterion and a charged membrane. We will call the 
interaction potential of such a system the Moreira-Netz 
potential. 

Our second system has for its interaction potential an 
attractive square well , which is a model for describing the 
adhesion of membranes by short-range attractive forces, 
studied e.g. in Ref. [31 . 

In our third system, the direct interaction is modeled 
by a Morse potential which has been used in studies of 
biological membrane systems and interactions of DNA 
molecules in solution mm [58] . The Morse potential 
consists of a repulsive exponential term which mimics the 
hydration and excluded volume repulsion of the mem¬ 
branes and an attractive exponential term (with twice 
the decay width of the repulsive term) which mimics the 
attractive long-raneg tail. 

Finally, we consider a system consisting of two mem¬ 
branes interacting via a generic soft potential consisting 
of a hydration term and a van der Waals term approxi¬ 
mated in the non-retarded planar limit [13 mg. 


1. Moreira-Netz potential 

Let us first consider the case of a pair of strongly 
charged membranes in solution with counterions of va¬ 
lence q in the intervening region, of charge sign opposite 
to the charge on the membranes. As the membranes 
are strongly charged, the system is in what is known as 
the strong coupling (SC) regime [5M5B] . where Poisson- 
Boltzmann (PB) theory and its fluctuation-corrected ver¬ 
sion break down. The SC regime is characterized by 
H > 1, where the coupling strength 5 = 2Trq 3 £^a s . Here, 
a s is the number density of charges on the surface of ei¬ 
ther membrane, and £ B = e 2 /47reeo^’BT is the Bjerrum 
length (e is the elementary charge). The Bjerrum length 
provides a measure of the interaction strength between 
a pair of unit charges at a given temperature T. The 
regime for which PB theory holds is known as the weak- 
coupling regime, and is characterized by S < 1. 

In the SC regime, the interaction between two equally 
charged hard planar surfaces is described by the Moreira- 
Netz free energy per unit area Mi IMl: 

/hard = 4.irk B T£ B cr 2 s - In , (3.24) 

where p = l/2nq£ B a s is the Gouy-Chapman length, 
which is inversely proportional to both the valence of the 


counterion species and the sur face c harge density of each 
plate. The first term in Eq. (3.24| reflects the electro¬ 


static attraction between the counterions and the charged 
plates, whilst the second term originates from the en- 
tropic pressure of the counterions confined between the 
two plates. The Moreira-Netz free energy was derived 
for a system where the charged plates are held fixed and 
are effectively rigid. 

We can generalize the Moreira-Netz free energy to the 
case of two thermally fluctuating membranes, where the 
membranes are allowed to equilibrate under an exter¬ 
nally applied osmotic pressure (thus we are considering 
a constant osmotic pressure ensemble). We replace i 
with £(xjJ, where the inter-membrane separation now 
depends on the transverse coordinate Xj_ = ( x,y ). As 
before, we write the separation as the sum of mean and 
fluctuating contributions: f(xj_) = £ 0 + 6£(x j_), where 
£o and 5£{x±) are thermal variables. The corresponding 
soft potential w is given by 


w = — 4 




= — AttIcbT/ liB&s 


ln£. 

2p 


2p 2p J 

6£(x±) (6£(x±)) 2 


2£ 2 


_£o_ _ (5l(xjJ \ 

2/i 2p ) 


(3.25) 


In the second step, we have Taylor expanded the log¬ 
arithm to quadratic order in 5£(x.±). On applying the 
Gaussian approximation of Eq. (3.12), we obtain 




(3.26) 


The variational free energy per unit area is then given by 
[cf. Eq. (pUgl)] 


/3/var = 


k B T 
128 no 2 + 

—4-7T p£ B a 2 


9c 2 a 2 3c 2 

+ Nia + 


K 


u 2 


, «o £o o- , 

n 2p 2/z 2f 2 J ^ ^ 


where P is the external osmotic pressure applied to keep 
the membranes at a constant average separation. Vary¬ 
ing / var with respect to a, we obtain 

-1/2 


= 1 + 


8 crJ, 


s^o 


9 c 2 q J 


k B T £ 2 
K 12 c 


(3.28) 


To form dimensionless variables, let us rescale a and £q 
in units of /to = qp, (since the Gouy-Chapman length p 
varies as q~ 1 ', we are interested in the effect of varying 
the valence q, so the basic lengthscale should be inde¬ 
pendent of q). Let us also define dimensionless variables 
t = k B T/ k, ( 0 = £q/po, 5/ = p 2 a s and? = cr/po, whence 
the above equation can be put in the form 


a 


2 


/ ~ \ -1/2 
Vi£l ( 8cr s lU 

12c y 9 c 2 q J 


(3.29) 
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p 


FIG. 3: Moreira-Netz potential. Behavior of £* = iq/pq 
(where l q is the saddle-point approximation to the equilib¬ 
rium inter-membrane separation and pq = 1/27t£b ct s ) as a 
function of rescaled osmotic pressure p = /3P /47r4o's [given 
by Eq. (3.32 l] and (inset) behavior of the rescaled rms fluc¬ 
tuation amplitude cf = o/pq as a function of £* [defined in 
Eq. (3.291] for t — 0.1, c = l and a s = l, studied for the fol¬ 
lowing three cases: (i) q = 2 (blue), (ii) q = 3 (green dashed); 
and (iii) q = 4 (red dot-dashed). For comparison, we have dis¬ 
played the behavior of the disjoining osmotic pressure (hor¬ 
izontal axis) due to counterions between two fixed charged 
plates as a function of the inter-plate separation (vertical axis) 
for counterion valences q = 2 (black dotted), q = 3 (cyan dot- 
dot-dashed) and q — A (orange dot-dashed-dashed) [obtained 
by differentiating Eq. \3.2A l], and also shown (cf. inset) the 
behavior of cf of a membrane interacting only sterically with 
a hard wall [black; cf. Eq. (2.261]. 


This formula holds strictly for the SC regime where S 
1. Substituting the value of a 2 into our expression for 
/var yields /soft, the effective interaction energy between 
two “soft” surfaces: 


where / = /3//4 Tip£ B <j 2 . In the saddle point approx¬ 
imation, the equilibrium separation £* is given by the 
solution to the equation df so ft(£q)/d£q = 0. We obtain 


P 


1 q 3 qc 2 

£* 2 + 8 a s (I*) 3 12c £* 

+ 1 

16 a s {£*) 3 


v ft 

V 1 + 9c 2 q 


(3.32) 


The last three terms on the RHS of Eq. (3.321 describe 
corrections to the equilibrium separation induced by the 
thermal fluctuations of the membrane. 


In Fig. § we display a plot of £* as a function of p for 
the case of thermally fluctuating membranes, where we 
have also plotted the case of fixed charged plates for com¬ 
parison. In the inset we show the behavior of a = cr/po 
as a function of £* for thermally fluctuating membranes. 
The behaviors are plotted for t = 0.1, c = 1, a s = 1, and 
three choices of the counterion valence: q = 2,3,4. As 
we see from the figure, in the SC regime a bound state al¬ 
ways forms at zero external osmotic pressure. Secondly, 
the bound state separation is smaller for larger counte¬ 
rion valences. Thirdly, at zero osmotic pressure the equi¬ 
librium bound state separation of thermally fluctuating 
membranes is consistently enhanced relative to the corre¬ 
sponding separation of the hard plates for the same coun¬ 
terion valences. Finally, at large values of £* the pressure 
curves tend to a saturation value of p = —q /2 for coun¬ 
terion valence q, identical to the case with no thermal 
fluctuations. Our prediction is thus that for multivalent 
counterion-mediated interactions, the equilibrium inter¬ 
membrane spacing for soft and rigid surfaces, everything 
else being the same, should differ and the larger the va¬ 
lency, the more they should differ. The predicted differ¬ 
ence in equilibrium spacing can be in excess of a factor 
of 2, see Fig. [3| 


2. Attractive square well potential 


/3/soft = 


3c lk B T / 8<j s £q \ 1/2 3c 2 


16 1 \ 


1 + 


9c 2 q ) 


+ ^ (3-30) 


Air p£ B o' 2 In -b (/3P + 2ir£ B a 2 )£ 0 . 

2 p 


The system still needs to be equilibrated with respect to 
£q. Let us define the dimensionl ess va riables: £q = £q/ pq 
and p = fiPjA- k£ b o 2 . Equation (3.30) then becomes 


/soft(4) = q(p + \)£o - In ^ + l 9 ?- 
^ A 16 a s £ 2 


+ 


3 qc\ft, 
32 a s £ 2 


/ 1 + 8 ° s £ 2 0 


9 c 2 q 


(3.31) 


Following Ref. [3Tj, we describe the binding potential 
by an attractive square well potential, and study the ef¬ 
fect that this square well potential has on the fluctuation 
and free energy behavior of a membrane using the varia¬ 
tional framework we developed in Sec. |III| As before we 
represent our system by a hard wall at z = 0 and a mem¬ 
brane whose surface is at a separation 2 = l/xjJ. The 
square well potential V(z ) is described by V(z) = V o for 
0 < z < b and V(z) = 0 for 2 > b. For an attractive 
potential Vq < 0. 

We represent the square well potential in terms of 
Heaviside functions 0, i.e., 

5?(*(x±)) = V 0 (O(£(x ± )) - 0(l( Xi ) - b)). (3.33) 
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In dimensionless form, the above equation becomes 


— — (fp-ir ' o 

£ 0 e 25 2 —(4-l)e 252 )+——- 


Vo\ 


t 

64 a 


,2^:3 


9c <7 


= 0. (3.38) 


The behavior of the effective interaction energy between 
the membranes as a function of £q can be determined 


from Eq. (3.36) subject to the contraint on a imposed 
by Eq. ( 3.38[ ). To determine the external osmotic pres¬ 
sure at which the membranes are maintained at a given 
separation, we differentiate Eq. (3.36) with respect to £ 0 


and set it to zero (which is equivalent to making the 
saddle-point approximation {£o)i ~ £q)- This yields 


FIG. 4: Attractive square well. Behavior of external osmotic 
pressure p (scaled in units of ksT/b) versus average separation 
l* for t = 0.1 and c, — 1, and the following four binding 
strengths: (i) Vo = 0 (black), (ii) vo = —0.04 (blue dashed), 
(iii) Vo = —0.06 (red dotted), and (iv) vo = —0.08 (green dot- 
dashed). Inset: behavior of dimensionless rms fluctuation 
<7 = cr/ b (where b is the well width) versus dimensionless 
separation £* = lo/b for the same values of t, c, and the 
above four binding strengths. 


The Gaussian approximation w a 2 is given by Eq. (3.12): 

dt o) 2 


W „2 = 


'-ooV^ira' 2 


V 0 (6(£) - &(£ - b)) e~ 


(3M) 


= Yi 

2 

The variational free energy [cf. Eq. (3.19)] is given by 


/var = y ( erf 


— erf 


+ 


(/cbT) 2 9/cbTc 2 ct 2 
128kct 2 + 8^ 


to-b 
y/2a 

Tc 2 


(3.35) 


+ Pto- 


We define a dimensionless separation distance i 0 = lo/b, 
dimensionless mean square fluctuation <7 = rr/6, dimen¬ 
sionless well depth v 0 = V 0 b 2 /V^ksT, dimensionless 
temperature t = k^T/n, and dimensionless external os¬ 
motic pressure p = Pb 3 /kBT. In terms of these quanti¬ 
ties, the dimensionless free energy is given by 


/ = /% 2 /var = (erf 


t 



— erf 



9c 2 ct 2 _ 

ns ^ + ^r + a +I ’ 


(3.36) 

C 0 Oi 0 

In the expression for the free energy above, a is not inde¬ 
pendent of io■ From Eq. (3.20), we find that a is related 
to to via 

V 0 




TUT* 


Up-6)2 

(£0 — b)e 2^2 — £ 0 e 


(/cbT) 2 9kBTc 2 a 


< 


= 0 


(3.37) 



For t = 0.1 and c = 1 the behavior_of the rescaled osmotic 
pressure p versus the separation £* and the behavior of 
the rescaled fluctuation amplitude d as a function of £* 
are plotted in Fig. [4] Note that the fluctuation ampli¬ 
tude is suppressed (relative to that of a purely steric sys¬ 
tem) for £q < b , whereas it is enhanced for £q > b. This is 
because for £q < b, the membrane effectively “sees” two 
repulsive (albeit one of which is finite) potential barriers 
and the repulsion has the effect of suppressing the ampli¬ 
tude of fluctuation, whereas for £q > 6, the membrane ef¬ 
fectively “sees” an attractive potential well which lessens 
the steric repulsion of the hard wall, and the fluctuation 
amplitude of the membrane is thus enhanced. For even 
larger values of £$, the fluctuation amplitude tends to¬ 
wards that of a purely steric membrane (represented by 
the black line), as it should. 

We note that the use of the Gaussian variational ap¬ 
proximation is reliable for sufficiently small well depths 
and/or membrane separations, a criterion of reliability 
being that the rms fluctuation behavior has to approach 
that of a purely steric membrane at sufficiently large 
membrane separations, as the square well potential is 
short-ranged. For larger well depths and/or larger sepa¬ 
rations, the application of the smooth Gaussian approx¬ 
imation to a sharp square well results in a certain os¬ 
cillatory behavior reminiscent of the Gibbs phenomenon 
(which can be observed in the rms fluctuation behavior at 
large membrane separations) and the predictions specif¬ 
ically for the order of the unbinding transition are not 
reliable. In fact the VGA-based formalism with attrac¬ 
tive square well potential and our steric potential predicts 
a discontinuous unbinding transition. This discrepancy 
with FRG-based approaches mi mi that predict a con¬ 
tinuous unbinding transition could arise from the rela¬ 
tive smallness of the co-dimension (which is unity) of the 
membrane [U] and/or from the fact that the hard steric 
potential is being approximated by the soft(er), long(er)- 
range interaction of a finite A, reflecting the size of a soft 
boundary region that depends on the microscopic details 
of the chemical make-up of the membrane m- 
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3. Morse potential 


Let us now consider the following Morse potential, 
which has been used to describe the condensation of DNA 
molecules in multivalent salts [571 [55] and to model the 
interactions of fluid membranes [7]: 

w(£) = w\{£) + zzq(^), (3.40a) 

u>i(£(x_l)) = Re- KD ^ X± \ (3.40b) 

w 2 (£(x± )) = Me.-5 SD<(xi ). (3.40c) 


The Morse potential has three fitting parameters: the 
strength of (short-range) repulsion R, strength of (longer- 
range) attraction A, and an inverse length-scale k b . Ap¬ 
plying the variational approximation to w a 2 [with the aid 
of Eqs. (3.12) and (3.40)], we obtain 


W a 2 = 


d£ 


— oc'/27 


,(Re~ KDe -Ae-* SD< ) e -1 ^ - 


= R e -^oto + k K O a2 (3.41) 


The variational free energy per unit area [defined in 
Eq. (3.19)] becomes 


/ var = Re-^o+k^ 2 _ J 4 e -)«D<o+l«^ 2 (3.42) 

, (faT) 2 9k B Tc 2 a 2 3fc B Tc 2 

+ 128kct 2 8£q + 8 £\ °' 


Defining dimensionless quantities u\ = [3k d 2 R, U 2 = 

Pk b 2 A, / = /3/«D 2 / var , £0 = kd4, ? = k d ct, p = Pk^P, 
and t = k B T / k 1 we can recast the above free energy in 
dimensionless form: 


/ 


zqe 


—^ 0 -t-^cr 2 _ 


u 2 e 2 ^+gff 2 


128a 2 


+ 


9c 2 a 2 

~W 


-p£ 0 . 


3c 2 

8£ 2 

(3.43) 


By varying / with respect to < 7 , we obtain a self- 
consistent relation: 


!(l e - 4 +j 5 2 _ -5A+5 52 _ ^ 1 ?£_ 

2 8 128a 4 8 £q 


0. (3.44) 


In the saddle-point approximation, the equilibrium sepa¬ 
ration (£ 0 ) is given by the solution £ 0 = £* to df /dig = 0: 


p = u 1 e~ e ° + ^ 2 


U2 ^ 0 + |5 2 
2 


Q r 2 /T 2 Q ^ 2 

+ (3.45) 

2^o 


In Fig. |5] the behavior of i* as a function of p and 
the behavior of a as a function of £* are plotted for 
c = 1, k = 10k B T , £* < yfS, and the following four 
cases: (i) zq = 112 = 0 (black), (ii) zq = 1,112 = 0.5 
(blue, dashed), (iii) zq = u 2 = 1 (red, dotted), and 
(iv) zq = l,zq = 2 (green, dot-dashed). Note that for 



FIG. 5: Morse potential. Behavior of l* = kd^o (where £q 
is the saddle-point approximation to (^o)i) as a function of 
external osmotic pressure p = pPnf) 3 an d (inset) behavior of 
rms fluctuation amplitude a = kb& as a function of separation 
£* for c = 1, k — IO&bT, £* <C \/S, and the following four 
cases: (i) Ui = 112 = 0 (black), (ii) ui = 1, 112 = 0.5 (blue, 
dashed), (iii) ui = 112 = 1 (red, dotted), and (iv) ui = 1,«2 = 
2 (green, dot-dashed). 


cases (ii) and (iii), the rms fluctuation of the membrane 
is smaller than the steric-only case (i) for £* < 1, as the 
repulsive interaction dominates in this range of separa¬ 
tions. At separations larger than £* ~ 1, the attractive U 2 
interaction dominates over the repulsive zq interaction, 
and the rms fluctuation is enhanced relative to the steric- 
only case, becoming larger for larger attraction strengths 
1 x 2 • At still larger separations the rms fluctuation of all 
four cases converge as the steric potential (which decays 
as (£*)~ 2 ) dominates over the exponentially decaying at¬ 
tractive tail. The rms fluctuation goes to zero at £* = 0 
owing to the hard wall constraint. For 112 > 0 the mem¬ 
branes are always bound at zero external osmotic pres¬ 
sure, and the bound state separation is smaller for larger 
values of zt 2 - 


4- Hydration and van der Waals forces 

Let us now consider the case of two uncharged mem¬ 
branes interacting via hydration and van der Waals 
forces, a generic case of interacting zwitterionic lipid 
membranes mm- For a pair of membranes of finite 
thickness 6, the interaction potential can be expressed 

as [171 HO] 


~/„n „ _W r 1 2 1 1 

w(£)- H e Ah l27r { e 2 (£+6) 2 + (£ + 25) 2 )' 

(3-46) 

where the first term describes repulsion due to hydration 
forces and the second term describes van der Waals at¬ 
traction of Hamaker strength W. For small membrane 
undulations, we can expand £(xj_) around £q to quadratic 
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FIG. 6: Hydration repulsion and van der Waals attraction. 
Behavior of dimensionless free energy / = dAii/var (with P = 
0 ) and (inset) dimensionless external osmotic pressure p with 
respect to rescaled separation £* = 4 /Ah for Ah = /MhAh = 
4.83, w = f3W/12n, 5 = <5/A H = 13.3, t = k B T/K = 0.0248, 
T = 270A', c = 0.255, and the following Hamaker strengths: 
(i) w = 0.04 (black), (ii) w = 0.04902 (orange dash-dash- 
dotted), (iii) w = 0.053 (red dot-dashed), (iv) ui = 0.061 
(blue disks), and (v) w = 0.068 (green diamonds). 


order in 5£(x j_). Using Eqs. (3.12 1 and (3.19), we obtain 
k B T 


Pfv ar = 


128 ^ + | + ^"^ + ^ 


9c 2 

w 0 


/3W 

47T 


1 2 , 1 

¥ 0 ~{^TW + W+W 

2 1 


_/3W 11 

”12^1^ (4 + <5) 2 1 (4 + 2<5) 2 


(3.47) 


Defining rescaled quantities t = k B T / k, 4 = 4/Ah, 8 = 
8/ Ah, Ah = p = /3PA^, w = PW/\2tt and 

/ = /3A^/ V ar, we can re-express the above equation in 
dimensionless form: 


/ = 


1285 2 8 £ 2 


3c 2 ~ _ 

-~r + A H e 


p 4 


+ • 


9c 2 


— 3w 


1 2 1 

3 (4 + 8) 4 (4 + 28) 4 

2 1 1 

w \ll (4 + 4 2 + (4 + 2?) 2 J 


(3.48) 


The relation between a 2 and 4 is given by Eq. (3.20), 
which yields 


A H e 


q i 1 ~2 

— «o+2 cr _ _ 


9c 2 


64<t 4 4^4 


(3.49) 


1 2 1 
—6 w < —-^^-b —- 

(4 + <5) 4 (4+ 2(5)4 


= 0. 


In Fig. [6j we study the behavior of the free energy and 
external osmotic pressure as functions of inter-membrane 


separation, for fixed hydration strength = 0.2 J m 2 , 
hydration lengthscale Ah = 0.3 nm, and bilayer thickness 
(5 = 4 nm. The critical point w o = w c at which dp/d£* = 
0 and d' 2 p/d(£*) 2 = 0 is determined numerically; we find 
that w c = 0.04902 for t = 0.0248 and c = 0.255. We 
also find the critical separation £ c = 26.9Ah and the 
critical rms fluctuation amplitude cr c = 4.76 Ah- For 
Hamaker strengths greater than w c the system exhibits 
phase co-existence similar to the one found in the case of 
the attractive square well (Sec. Ill A 2). The threshold 
Hamaker strength at which the membranes undergo a 
discontinuous and complete unbinding transition at zero 
external osmotic pressure is estimated to be Wd = 0.061 
(the blue disks in Fig. [6]). 

Whereas FRG-based approaches (20125 predicted a 
continuous unbinding transition for a pair of steric mem¬ 
branes interacting via van der Waals (vdw) and hydration 
forces, our VGA-based formalism with our steric poten¬ 
tial predicts a discontinuous unbinding transition. Just 
as we concluded for the short-ranged square-well poten¬ 
tial in this case too the VGA-based formalism as imple¬ 
mented here does not reliably predict the order of the un¬ 
binding transition when compared with the FRG-based 
approaches 120125 for the same reasons as already in¬ 
voked above. 


IV. SUMMARY AND DISCUSSION 


We have proposed a self-consistent theory for studying 
the interaction between a pair of mutually impenetrable 
and thermally undulating fluid membranes giving in the 
process a definitive and consistent form to the previous 
partially successful attempts in the same direction 0 ®. 

We have implemented the steric constraint via the 
Panyukov-Rabin representation of the Heaviside func¬ 
tion. For a pair of membranes of bending stiffness k and 
cross-sectional area S interacting exclusively via steric 
and fluctuation forces, and separated by a mean distance 
4 , we have derived a closed-form expression for the steric 
potential per unit area V s [see Eq. (2.29)]: 




3 k B Tc 2 k B T r 9k B Tc 2 S 2 \ 

8£q 87 tS \ 4k£q J 


, k B T j 9k B Tc 2 f ^ .4 V2(4 K y/ 4 £ 0 ) 

+ 4tt Y 4n£\ I an V (9k B Tc 2 y/ 4 VS ) 


+ tan 


/ V2(4k)V% 

V (9k B Tc 2 y/ 4 Vs J) 


(4.1) 


This has two contributions: one that is induced by zero 
mode fluctuations of the membranes and one that is in¬ 
duced by thermal bending fluctuations. At small sep¬ 
arations 4 •C y/S, the bending fluctuation-dependent 
part scales as £// 2 , and crosses over to £g 4 scaling for 
large separations 4 ^ VS. On the other hand, the zero 
mode-dependent part of the steric potential always scales 
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as£ 0 1 2 3 4 5 6 7 8 . Concomitantly we also derived a closed-form for¬ 
mula [see Eq. (2.23)] for the rms undulation amplitude 


f 2 

a 2 = AIL 
12 c 


k B T 


1 -tan 

7T 


K 


3cS V k B T 


(4.2) 


which has the following asymptotic behavior: for £q -C 
y/S, a scales linearly with £ 0 , whereas for i 0 3> V~S, <r 
saturates at the order of k B TS/n. The rms fluctuation 
amplitude becomes small at low temperatures T and/or 
large bending stiffnesses k. We believe that our result 
refines and substantiates the Ansatz er 2 = jll 2 first pos¬ 
tulated (on the basis of heuristic arguments) in Ref. [2]. 

To investigate fluid membrane systems that experi¬ 
ence interactions of non-Gaussian form, we have adapted 
the Feynman-Kleinert version of the variational Gaus¬ 
sian approximation (VGA) to the case of fluid mem¬ 
branes subject to our effective steric potential V s . We 
have applied this VGA approach to four different types 
of potential: (i) the Moreira-Netz potential for a pair of 
strongly charged membranes with an intervening solu¬ 
tion of multivalent counterions, (ii) an attractive square 
well, (iii) the Morse potential, and (iv) a combination of 
hydration and van der Waals interactions. 

In the first case we make a prediction that, every¬ 
thing else being the same, the multivalent counterion- 
mediated interaction measured between hard and soft 
surfaces should display a substantial difference in the 
equilibrium spacing, with soft surfaces displaying larger 
equilibrium spacing. This difference should increase with 
the valency of the counterions and could easily reach a 


factor of 2 for high valency counterions and should thus 
be eminently measurable in the planned osmotic-stress 
experiments with lipid membranes in the presence of mul¬ 
tivalent salts [55]. We furthermore note here that our 
results for multivalent salts imply also a pronounced ef¬ 
fect on the estimated values of the bending rigidities that 
could be extracted from the linewidth of the X-ray scat¬ 
tering intensity in this type of experiments, irrespective 
of whether they directly renormalize the bending rigidity 
of lipid membranes or not |60H62j . 

Within the same VGA approach we also analyzed 
in detail the three other cases of the coupling between 
long-range interactions and conformational fluctuations, 
showing the versatility and usefulness of our approach in 
the context of widely differing types of long-range inter¬ 
actions in the regime of intermediate separations, which 
are also the typical length scales of biologically relevant 
systems. 

As future projects in the same general direction, our 
self-consistent theory can be employed to study multi- 
component membranes, membranes with non-zero sur¬ 
face tension, and tethered and polymerized membranes. 
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